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I , Low energy phenomena involving two nucleons can be successfully described using 

^ 1 effective field theory. Because of the relatively large expansion parameter, it is only 

3 ' at next-to-next-to-leading order (NNLO) where one can expect to see agreement 

with experiment at the few percent level. The first part of this talk will focus on ra- 
ff^ ' diation pion effects, which first appear at NNLO. The power counting for radiation 

pions is simple for center of mass momentum p ~ \/~Mrn^ = Qr, the threshold for 
pion production. We explain how graphs calculated with the Qr power counting 
^ ' scale for p ~ m,r • The Q^ radiation pion contributions to nucleon-nucleon scatter- 

^~v ' ing are suppressed by inverse powers of the S-wave scattering lengths. However, we 

point out that order Qf radiation contributions might give a NNLO contribution 
for p ~ m,r . In the second part of the talk, results for the potential pion and con- 
tact interaction part of the NNLO ^Sq phase shift are presented. We emphasize 
the importance of eliminating spurious poles in the expression for the amplitude 
at each order in the perturbative expansion. Doing this leaves a total of three free 
parameters at NNLO. We obtain a good fit to the ^ So phase shift. 



Introduction 



This talk focuses on higher order calculations in the low energy effective field 
theory for two-nucleon systems. In particular, we will be discussing nucleon- 
nucleon scattering at next-to- next-to-leading order (NNLOl in the expansion 
recently proposed by Kaplan, Savage and Wise (KSW)El'H Observables are 
expanded in powers of Q/A, where Q is either p, the three-momentum of the 
^ I two nucleons in the center of mass frame, or m^^. A is the range of the effective 

theory. The nuclear S-wave scattering lengths (denoted by a) are very large 
so that powers of pa have to be summed to all orders at each order of the 
Q expansion. This requires a novel power counting in which the leading 4- 
nucleon operator with no derivatives is treated nonperturbatively To make 
this power counting manifest in dirnensional jregularization it is necessary to 
use subtraction schemes such as PDSEla or OSO'cl, but predictions of the theory 
are manifestly scheme and scale independent order by order in Q. Higher 
derivative operators and pion exchange are treated perturbatively 

Various estimates of the range of the theory exist. One estimate comes 
from examining pion exchange ladder graphs, where each additional loop gives 
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a contribution of order p x Mg'\/{8TTf'^) = p/K^N, where Ajvat — 300 MeV. 
Another possibihty is that rup or the threshold for A production sets the scale 
for the breakdown of the effective theory, implying a range for S-wave scattering 
~ 700 MeV. It has been suggested ErQ that two pion exchange contributions 
to the nucleon-nucleon potential may be become important for momenta of 
order 400 MeV. These considerations point to a range somewhere between 
300 MeV and 700 MeV. Therefore, for p ^ 111^^, the expansion parameter, 
Q/A, is between 1/2 and 1/5. Because the expansion parameter of the theory 
is rather large, low order calculations in the effective theory cannot be expected 
to reproduce phase shift data as accurately as potential models with many 
parameters. 

Many observables have been coninuted to NLO in the KSW expansion, in- 
cluding nucleon-nucleon phase shiftsPJiO Coulomb correctioas to proton-proton 
scatteringP electromagnetic forna. factors ot-jche deuteronp deuteron polaria* 
abilitiesp^ pratpn-proton fusionpJ np — > d'j^ Compton deuteron scatteringpJ 
and vd -^ i>dt3 Somejaf these calculations are reviewed in the talk by Martin 
Savage in this volumellj One typically finds errors of order 30 — 40% at leading 
order and 10% at NLO. This is consistent with Q/A « 1/3, or A « 400 MeV. 
This suggests that at NNLO, effective field theory calculations of low energy 
processes in the two body sector should agree with data at the few percent 
level, approaching an accuracy comparable to that of potential models. It is 
for this reason that extending calculations to this order is an important part 
of the effective field theory program. 

In the first half of this talk, we will discuss radiation pion effects, which 
first appear at NNLO in calculations of nucleon-nucleon scattering. The power 
counting of KSW has to be modifiedjJxi the presence of pion radiation because 
a new scale, Qr = y/Mm^, appearsllj For power counting radiation pions, it 
is simpler to take p ~ Qr and count powers of Qr rather than Q. We give 
a procedure for determining how a Q" correction scales with Q for p ^ m.,^. 
The order Qj? radiation pion contribution to nucleon-nucleon scattering turns 
out to be suppressed by powers of 1/a. This is actually a consequence of 
the invariance of the leading order theory under Wigner's S'[/(4) spin-isospin 
symmetryEj Winner symmetry is discussed in detail in the talk by Mark Wise 
in this volumeLfl The order Q^ radiation pion corrections can give an order Q 
contribution. In the second half of the talk, we present results of a NNLO cal- 
culation of nucleon-nucleon scattering in the ^Sq channel. It is emphasized that 
coupling constants of the theory must be treated in a Q expansion. The pa- 
rameter space is constrained by the requirement that perturbative corrections 
do not shift the location of the pole in the amplitude and by the solutions of 
renormalization group equations. Once these constraints are imposed, a three 



parameter fit to the ^5*0 phase shift is demonstrated which has < 2% accuracy 
at p ~ rra^ and also reproduces the data weU for higlpr—mornenta. A similar 
calculation is discussed in the talk by Gautum RupaktSa. 

Radiation and Soft Pions 

The Lagrangian for the theory of nucleons and pions is 

£, = ^Tr (a^E a^St) + ^ Tr(m,I] + m,St) + '-^ N^a.i^d,^ - ^d^^N 

+ N^ (*^o + ^) ^ - Q{N^PtNyiN^P,^N) (1) 

+ % {N'^ P^' N)\N^ P^' ^^ N) + h.c. 
8 

- D^ujTT{m^){N^Pi'N)''{N^PlN) + ... , 

where operators relevant at NLO are included (and isospin violation is ne- 
glected). Here gA = 1-25 is the nucleon axial- vector coupling, S = ^^ is the 
exponential of pion fields, / = 131 MeV is the pion decay constant, m^ = 
(^"^9^ + C^'^qC^)/^, where niq = diag(?n„, tti^) is the quark mass matrix, and 
TO^ = w{mu + rnd). The matrices P/ project onto states of definite spin and 
isospin, and the superscript s denotes the partial wave amplitude mediated 
by the operator. This talk will be concerned only with S-wave scattering, so 
s = S (for ^5*0) or T (for ^Si). This notation will be omitted when it is not 
necessary to distinguish between the two channels. 

In the KSW power counting, the Cq operator is treated nonperturbatively 
and graphs with a single pion exchange (dressed with Cq bubbles) first appear 
at NLO. Loop graphs with pions contain three different kinds of contributions, 
which are called potential, radiation and soft. The three kinds of pion are 
characterized by different energy (go) and momentum ((f): 

potential qo ^ q /M 

radiation qo ^ \q\ ^ m^ 



soft go ^ \q\ ^ Qr ^ y/Mm^. 

As stated earlier, when calculating radiation and soft contributions we take 
p ~ Qr rather than p ^ Q. The three contributions will differ in size and it is 
necessary to devise a power counting which correctly takes this into account. 
Before giving the power counting we will illustrate how these contributions 
arise with a few illustrative examples. 



Figure 1; One loop graph with a Co and pion that has both potential and radiation contri- 
butions 

Consider the one loop graph shown m Fig. n, whose contribution to the 
ampHtude is proportional to 

g\Co rA 1 1 f 



When the qo integral is performed via contour integration, the integral receives 
contributions from both pion and nucleon poles. If a nucleon pole is taken then 
\qo\ = E/2 - {q -pY/2M = {2q ■ p - q'^)/2M, where in the last step we have 
used E — p /M. Since |p|,|(f| ^ M, go ^ \q\, we can expand the pion 
propagator: 

1 _ 1 ll , 



qi — q — rnl q + m^ [q + mIY 

A pion is referred to as a potential pion whenever the energy dependent piece 
of its propagator is treated perturbatively. In the KSW power counting, |p | ~ 
TOtt ~ M_R '^ \q\ ^ Q- Each nucleon propagator gives a factor of M/Q^ 
since E ^ q^ ^ Q"^ jM. The measure di^q ~ Q^ /M ^ and in a scheme with 
manifest power counting, Co ~ l/{M^ji) ^ 1/{MQ). Using this counting it is 
straightforward to see that this graph is order Q^, i.e., it is a NLO contribution. 
The first correction from the expansion in q^/{q + mj) is suppressed relative 
to the leading potential contribution by Q^/M^, and so is N^LO. 

There is also a contribution from the pion pole. In this case q^ = q -\- 
m^. In the nucleon propagators, the factors of {2q ■ p — q )/2M <C qo and 
must be treated perturbatively. With the KSW power counting, the nucleon 
propagators in the graph in Fig. 1 are ±l/go '^ 1/Qj a-nd the loop measure 
scales as Q**, so the graph is ~ Q. Therefore, this radiation pion contribution 
first appears at NNLO. 

While KSW power counting works for the graph in Fig. |^, it fails for other 
graphs with radiation pions. As an example, consider the graph in Fig. ^ which 




Figure 2: A radiation pion graph with n internal Co bubbles. 



contains n nucleon bubbles inside a radiation pion loop. The loop integral 
in this graph vanishes if the pion pole is not taken so there is no potential 
pion contribution. Emission of the radiation pion in these graphs changes the 
spin/isospin of the nucleon pair. Therefore, if the external nucleons are in 
a spin-triplet (singlet) state, then the coefficients appearing in the internal 
bubble sum are Cq (CJ). For definiteness, consider nucleon- nucleon scattering 
in the ^5*0 channel. The contribution from the graph in Fig. is: 
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The go integral is closed around the one pion pole above the real axis so go ^ 
\q\ ^ Q- In the fco integrals a nucleon pole must be taken. In the KSW 
power counting, fco ^ Q^/-^, \k\ ^ Q- The graph then scales as ((5/M)"+^. 
This suggests that graphs with nucleon bubbles inside the radiation loop are 
suppressed relative to the one loop radiation pion graph. However, explicitly 
performing the qo and fco integrals gives 
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The size of the loop momenta fc in the nucleon bubbles is 
p < ^Mm-n- The integral will be dominated hy q ^ m^ 



~^ \/Mm^ even for 
so the graph will 
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(Recall A^ = iirf). For /i^ r^ Q r^ rriTr, we see that graphs with additional 
bubbles are enhanced, contrary to the KSW power counting. The sum over 
graphs with an arbitrary number of bubbles is /i/j independent and the correct 
estimate for the size is obtained when hr ~ \fMrn^. At this scale, these 
graphs and their sum are of order •mj'^ l{M'^l'^ls2^ = (5^/(M^A^). 

In the KSW power counting, one assumes the loop momentum in the nu- 
cleon bubbles is dominated by /c ~ Q, fco ~ Q^ /M . However, when the nucleon 
bubbles are inside a radiation pion loop, the energy flowing into the nucleon 
bubbles is actually order TTi^r, and therefore kP ~ m-^ and k ~ \/Mra^ = Qr- 
This scale corresponds to the threshold for on-shell pion production. 

In general, radiation pion graphs will depend on p, TOtt and M in a compli- 
cated way, making them difficult to power count if one takes p ^ rriT^. However, 
the natural scale for loop momenta with radiation pions is Qr, and the power 
counting simplifies considerably if we consider nucleons scattering at p ~ Qr- 
Later we will discuss what happens as p is lowered back down to ttt,^. 

Power counting at the scale Q^ is straightforward. We take p ~ /i/j ~ Qr- 
The C2n scale with /i/? exactly the same way as in the KSW power counting, 
C2nP^^ ^ P^"/Mfl^^ ^ Qr^^- A radiation pion propagator gives M'^/Qf, the 
pion nucleon coupling gives Ql/Al . Nucleon propagators scale like M /Q^. In 



a radiation loop go 



so the loop measure d q ~ Q^/M . The 



measure of a potential loop scales as Q^/M. Using this power counting it 
is straightforward to show that all graphs with one radiation pion and an 
arbitrary number of Co's scale as Q'^/{M^KV). These graphs are shown in 

Fig. I 

It is interestirLg^to examine the result of evaluating some of the diagrams 
in Fig. I explicitly^ {fi^ = ^^^e"^^): 
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Figure 3: Leading order radiation pion graphs for A'^A'^ scattering. The solid lines are nucle- 
ons, the wavy lines are radiation pions and SM, 5Z are the mass and field renormalization 
counterterms. The filled dot denotes the Co{ij.r) bubble chain. There is a furtjier field 
renormalization contribution that is not shown, but is included in the calculationr^ 



Ji and I2 are hypergeometric functions. The 1/e poles are cancelled by inser- 
tions of a D2ml counterterni. The leading order amplitude 
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scales as '- 


-^ l/{AIp), so we see that Eq. (2) has terms proportional to 




/m,x3/2 ml m^/' 




Km) ' Mp M^/^p^ ■ 



For p '^ Qr these terms scale as Q^/M^, as anticipated by the power count- 
ing. At p ~ TO^ ~ Q, these terms scale like {Q/Mf/^,Q/M, and {Q/MY^^ 
respectively. Bubble sums which do not appear inside radiation loops will be 
referred to as external bubble sums. These bubble sums are responsible for 
the factors of p in the denominators, and thus the enhancement of some terms 
at low momentum. Graphs with two external bubble sums have terms that 
are enhanced by Qr/Q^ ^ I/'tItt (and Qr/Q ^ I/^/^Itt), while graphs with 
one external bubble sum have terms enhanced by l/y/m^. Individual graphs 
like b) have parts that scale differently with Q. Terms which scale like Q^^^ at 
low p are actually larger than NNLO in the Q counting. The Q^^^ contribu- 
tions come from graphs 6),e) and /), and cancel when the graphs are added 
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together. Presently it is not known whether this canceUation occurs for some 
reason or is merely an accident. 



The sum of all Qf. graphs is 



-\4Trf)^ 4tt \aS qT 



(4) 



In 



where the ln(/i) dependence is cancelled by D^^^i). (For the '^S\ channel, 
the result is the same as above with as <-> ot-) The final answer turns out 
to be much smaller than anticipated by the power counting. For p '^ Q^, 
the first term is suppressed by a factor of XjQri^ja^ — 1/a^)^ the second 
by l/Q'^{l/a^ — 1/a^Y . This suppression occurs because the radiation pions 
couple to a charge of Wigner's S'C/(4), which is a symmetry of the leading order 
Lagrangian in the limit 05, a^ — > cx) (or as = aT)tZI'E3 The order Q^ radiation 
pion graphs are a tiny correction to the S-wave scattering amplitude. 

The next important radiation pion contribution comes from graphs with 
one insertion of a C2P^, I?2"t-^, or G2 operator □ or a potential pion, and 
one radiation pion with an arbitrary number of Cq's. Power counting these 
graphs gives Qf/{AI^A'^A), i.e. these are suppressed by Qr/A relative to the 
leading radiation pion graphs in Fig. pi Note that Qr = 360 MeV, so for the 
most pessimistic estimates of A, the Qr/ A expansion does not converge. If 
this is the case then the radiation pion contribution is incalculable. This is 
true of radiation contributions even when we scale down to p ^ tjIt^. How- 
ever, at the low momenta where the theory would be applicable, the radia- 
tion pions could be integrated out. For example, the hypergeometric function 
/2(-E/to^) = /2(p^/(Mr?T,7r)) could be expanded in a series in p^ and the effect 
of each term absorbecl|^ into the definition of a C2n- Assuming the radiation 
pion contribution is computable, the Qf correction is almost certainly larger 
than the Q^ correction. Since the C2P^, D2m^ operators and potential pion ex- 
change do not respect Wigner's SU{4), there will be no suppression by factors 

of l/(aQ.)- 

An important issue which needs to be addressed is how the radiation pions 

graphs scale as p is lowered from Qr to m^. We saw that the Q^ graphs had 
pieces that scaled as Q^^^, Q, Q^^^, ■ ■ ., for p ^ m^r, and that this enhancement 
can be understood by counting the number of external bubble sums. In order 
to know which radiation pion graphs to include at a given order in the KSW 



''Unfortunately, the resulting theory below the scale Qr would no longer respect chiral sym- 
metry. Operators involving a different number of pion fields could have different coefficients. 



power counting, we must know how a Q" correction scales with Q for p ^ 111^,. 

It turns out that an order Q" calculation is sufficient to determine the order 
gn/2-1 j.gg^i^ 

To see this first consider the Q expansion of p cot (5: 
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pcot i5 is real and an analytic function of p^ near p = 0. This will be true order 
by order in Q so: 

-72~ = /" ^ ^0 = hA_i , 



Al At, 



All -4-1 
A2 2A0A1 A^ 



yn. 1 J^ 1 -rX -1 



= /i ^ ^1 = hAli + flAU , 

= h ^ A2 = f2A\ + 2,fo.fiAU + f^Ati , 



where the /„ are real functions of p which are analytic about p'^ = 0. We see 
that the general form of a higher order amplitude is powers of A-i multiplied 
by functions of p. The crucial point is that the function multiplying the Ali 
is the only new contribution. The coefficient of A"Li, n > 2, is determined by 
lower order amplitudes. In the Q expansion of p cot (5 the latter contributions 
will cancel. 

This generalizes to the Qr expansion of radiation pion graphs, the only 
difference being that the radiation pion contribution starts out at Q^, while 
the potential pion starts out at Qj?. A Q" radiation pion correction to the 
amplitude will be of the form: 

An = A_ifn,2 + A'_i,fn,3 + . • . + A_i ,fn,n-l- 

Again, the fn,m are analytic about p^ = and all the fn.m except for /„^2 
will be determined from lower order amplitudes. Since An ^ Qr and A-i ^ 
l/{Mp), fn,m ~ Qr^"^ for p ~ Qr- To Understand how fn,m scales with Q as 
p is lowered to m^r, note that without loss of generality, fn,m can be written as 






y*(lXJ*x y*CXXiJ*x x*U*(i)*x 



Figure 4: Example of order Qf. graphs that have three external bubble sums. 



where the elhpses denote momentum dependence that involves scales other 



than Qr, and A = A^., A, or M. For p 



the ellipse denote dependence 



on the dimensionless variables p/niTr, pa, and p/A. For p ~ TTijr, p/y/Mm^ 
{Q/My'^ and the function fn,m can be expanded in its first argument: 






A^A 



01^ 

M 



1/2' 



The leading term scales as '--^ Qn/2-m/2 jt^j. ^ ^ ^^^ ^ q ^{j^^q only the ?ti = 2 
term actually ends up contributing to pcot 5, the new contribution at Q" scales 
like Q"/^^^ (plus subleading terms) for p ~ m^r. This is consistent with the 
result of the Q^ calculation, where the largest contributions from individual 
graphs scaled as Q^^^ . A cancellation between graphs resulted in this contribu- 
tion vanishing. The remaining terms scale as Q, Q"^'^, . . . (counting 1/a ^ Q). 
The Q" radiation pion contribution could have a Q^'^ contribution from the 
A^ fn.n-i term. But this contribution is determined by the Q'}, amplitude 
which vanished; so there will be no Q^'^ contribution from any radiation pion 
graph. For example, at Q^ the terms proportional to A^i comes from graphs 
such as those shown in Fig. ^ These graphs factorize into two pieces which 
are lower order and it is easy to see that the same cancellation between the 
Q^/^ pieces of the graphs b), e) and /) will also occur at Qf. 

Since Q"/^^^ = Q for 7i ~ 4, the Qf. radiation pion graphs may have a 
contribution that is NNLO for p ~ rn^. We have checked that graphs with 
one insertion of the D2m\ operator do not give rise to such a contribution, 
but a calculation of the remaining Q^ graphs has not been performed. These 
graphs need to be computed in order to obtain the complete NNLO amplitude. 
Unfortunately, graphs with one potential and one radiation pion are numerous. 
Higher Q" amplitudes may have a contribution which scales as Q for p ^ ?tt,^ 
from the yl"7 fn.n-2 term. But these will cancel in the Q expansion of p cot (5. 
Note that a calculation of the order Q^ graphs would be necessary to determine 
the order Q^^'^ terms. 
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a) b) c) d) 

► N S ► 



s: X 



Figure 5: Examples of one-loop graphs which have soft pion contributions. Graphs a)-d) 
also have a radiation pion contribution, while in addition graph a) has a potential pion 
contribution. 



Finally, we briefly discuss soft contributions. These first arise in box- type 
diagrams with two pions, such as those shown in Fig. ^. Dressing these graphs 
on the outside with Cq bubbles gives further diagrams of the same order. Un- 
like radiation graphs, which are dominated by loop momenta \q\ ~ m^, the 
box-like diagrams receive nonvanishing contributions from \q\ ^ Qr- It is 
this latter type of contribution which is called soft. In soft loops the measure 
df^q ~ Qf instead of Qf/M^. For nucleon propagators in soft loops, the loop 
energy is always greater than the nucleon's kinetic energy, soJche propagator 
is static (like heavy quark effective theory propagators, seeEj orEj). These 
propagators scale as 1/Qr- Pion propagators scale as l/Q^ and pion-nucleon 
vertices give a factor of Qr- Power counting the graphs in Fig. || gives ~ Q^ 
for the soft contribution and ~ Q^/M'^ for the radiation contribution. The 
contribution frorp-fjach .regime can be separated .using the method of asymp- 
totic expansions E^^ESEj and explicit evaluation t3 of the graphs verifies the 
power counting appropriate for each type of contribution. At p ^ Qr, the soft 
contribution actually dominates all radiation pion graphs we have considered 
so far. However, for p ~ TTi^r these soft graphs are order Q^, and therefore are 
not enhanced by the scaling down to p ^ wLtj . The leading order soft graphs 
are N^LO in the KSW power counting. 

NNLO Calculation of the ^^o Phase Shiftg 

In this section, we present a partial NNLO calculation of the^5o phase shift. 
A more detailed analysis will be given in a future publicationcj The first piece 
of the Q^ radiation pion contribution in Eq. (0) is order Q. The second term 
is order Q^/'^ and is not included in the NNLO calculation. The relativistic 
corrections are computed in RefED and shown to be negligible. This is because 
they are suppressed relative to the leading order amplitude by (Q/M)^ rather 



'^The work presented in this section was done in collaboration with Sean Fleming. 
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than (Q/A)'^ and so are smaller than other NNLO corrections. The calculation 
here includes order Q contact interactions and potential pion graphs. The 
NNLO calculation is incomplete because of the omission of the Q^ radiation 
pion graphs, as discussed in the previous section. The pieces of the NNLO 
amplitude which are included are expected to scale as Q/{Mh.^) while the Q^ 
radiation pion graphs are expected to scale as Q/{K^K) for p ~ m^. Since 
A < A^ w A/, the contribution from Q^ radiation pion graphs may be smaller 
than what has been included. 

Since we are only interested in the ^Sq channel, only s — S operators are 
relevant and this superscript will be omitted in the following discussion. At 
NNLO, the following terms are added to the Lagrangian in Eq. ™): 



64 



{N^P,N)\N^P,V N) + h.c. + 2{N^ P,V N)\N'^P,y N) 



8 



w2 

{N'^P,N)^N'^P,\/ N) + h.c. 



Da 2 



|Tr2(TO«) + 2Tr[{m^f]\{N^P^N)\N^P^N) , (4) 



where only terms relevant for the phase shift are included and isospin violation 
is neglected. All calculations presented in this section will be done in the PDS 
renormalization schemeo'EI with spin and isospin traces done in four dimensions. 
There are six coefficients that appear at NNLO: Cq, which is present in the 
leading order calculation, C2 and D2, which first appear at NLO, and C^E4^, 
and Di. It is important that the coupling constants are expanded in Q: H 

Co — > Co + Co,0 + Co,! 

C2 — > C2 + C2,-l 

D2^ D2 + D2,-i . (5) 

The first piece of Co is treated nonperturbatively (i.e. Co ~ Q^^), while 
Co.o ^ Q'^7Co,i ~ Q. Solving the renormalization group equations (RGE) 
for the couplings perturbatively ensures that the amplitude is /i independent 
order by order in the expansion. Therefore, theoretical expressions for physical 
quantities, such as the scattering length or the location of the pole, are always 
/i independent. Physically, the Q expansion of the couplings is a consequence 
of the fact that higher order loop graphs with pions can renormalize the short 
distance operators at different orders in Q, and therefore different values of the 
couplings will be obtained at different orders in the expansion. 

When a coupling is expanded in Q and its RGE solved perturbatively, a 
new constant of integration is obtained for each term in the expansion. For 
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example, the RGE for Cq is: 



u^ r - ^^^ (r + ^^ 



(6) 



After the perturbative expansion of Co this becomes 
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An 
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-M + 7 
1 
'-C-0 



Q 



0,0 



Co,o 



An ^° 2/2 



2/2 



Mk' 

An 



C' 



(7) 



Three constants of integration appear: 7, k and k'; one for each term in the Q 
expansion of Co- For consistency we must assign these constants a Q counting. 
For instance, Co,o has a term in its solution kCq ~ k/^'^. Since Co,o ~ Q^, 
K must be of order Q^ . This reflects the fact that k, is intrinsically small. In 
the theory without pions, k = 7 — 1/a. As we will see below, the values of 
7, K, k' may be fixed by demanding that the amplitude has the correct pole 
structure. In this case, k « ro/(2a2) ~ Q'^. The results in Eq. (0) can also be 
obtained by solving Eq. (g) exactly and expanding the result (including the 
constant of integration) in Q. Because of the perturbative expansion of the 
couplings in Eq. (g|) there are ten constants of integration at NNLO. However, 
the NNLO amplitude will depend only on six independent linear combinations 
of these constants. There are two further constraints on the number of free 
parameters: 1) at this order, C4, -E4 and 1^4 are determined entirely in terms 
of lower order couplings; 2) spurious double and triple poles in the NLO and 
NNLO amplitudes must be cancelled in order to obtain a good fit. 

The fact that C4 , -E4 and D4 are determined in terms of lower order cou- 
plings is a consequence of solving the RGB's and applying the KSW power 
counting. For instancecl, the RGE for C4 is: 



^|l^^ = ^('^°^^ + ^^') 
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which has the solution 

where p is a constant of integration. In the theory without pions, p is pro- 
portional to the shape parameter, which is '-^ Q° in the KSW power counting. 
It is also reasonable to consider p ~ Q^ in the theory with pions, since p is a 
constant of integration in the RGE for the lowest order term in the Q expan- 
sion of d- Therefore, C|/Co ~ Q-^ while pC^M/(47r) ~ Q-^. The second 
term is subleading in the Q expansion, and should be omitted at NNLO, so 
C4 = CI/Cq. Similar relations for E^, D4 hold at NNLO: 

2C2D2 ^ D^ 

It is interesting to note that these relations arise even though there are ln(/i^) 
terms in the amplitude which contribute to the beta functions for E^ and D^. 
Because of the nonperturbative treatment of Co, spurious poles arise at 
higher orders in the expansion. The leading order amplitude A-i has a simple 
pole dXp — ij. The NLO calculation is proportional to A'^i, and therefore has 
a double pole, while the NNLO amplitude has terms proportional to A^i and 
A'Li- To obtain a good fit at low momentum, parameters need to be fixed so 
that the amplitude has only a simple pole at each order in the expansion. This 
requires that A-i have its pole in the correct location and that the residues 
of the spurious double and trjple poles vanish. This requirement leads to the 
following good fit conditions: 



A-i 



0, ^' 



p=p' 



AU 



= 0, 



Ax 



AU 



= 0, (8) 

p=j}' 



where p* is the location of the physical pole. The second condition first appears 
at NLO, the third at NNLO. The residue of the triple pole in Ai is cancelled 
by the second equation in Eq. (||). The first equation results in 7 = —ip*, 
while the other equations give constraints which eliminate two of the remaining 
parameters. Eq. (ra) will also apply in the ^Si channel. 

The graphs evaluated in the NNLO calculation are shown in Fig. 6. The 
final result is surprisingly simple. The amplitude up to NNLO is: 

An 1 

IVl J + tp 
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Figure 6: Order Q contact interaction and potential pion graphs for the So channel. The 
shaded circle denotes a Co bubble sum. At this order the first six graphs cancel each other 
as explained in the text. 
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Using this amplitude it is easy to verify that the S-matrix is unitary to the 
order we are working. The six Uncarly independent constants appearing in the 
amphtude are 7, Ci, C2) Csi Ci^ Cs- By definition Ci — C5 are dimensionless. They 
are given in terms of couphng constants in the Appendix. 

For the ^5*0 channel, the location of the pole is determined by solving 



a 2 



ip 



0. 



(10) 



Adding the shape parameter correction to Eq. (|lO| ) changes the location of the 
pole by less than 0.01%, so this approximation for the physical pole location 
is sufficiently accurate. Eq. (|lO| ) has a solution for p* — — i 7.877 MeV, which 
fixes the single LO parameter, 7 — — 7.877 MeV. The NLO good fit condition 
relates the constants Ci and (2, 



Ci7 + (2^2 + - ^^ 



=^C2 = -tC 



Mg\M 



^log 



47r 87r/2 



(- 


472 
ml 


27 


27^ 


m-TT 


ml 



- tanh ^ —^ 



27 
m-, 



= 

(11) 



We can use this equation to eliminate C2 in favor of Ci , leaving one new parame- 
ter in the fit at NLO. This good fit condition gives non-trivial m-^ dependei 
the perturbative contributions to Cg (such as k) as emphasized in Refs.E 
Ci is fixed by doing a weighted least squares fit to low momentum data. Note 
that Ci appears in the good fit condition multiplied by 7^. Therefore, the 
value of C,2 is insensitive to the value of ^ obtained from the fitting procedure. 
To a good degree of accuracy we can ignore Ci in Eq. (ll|), and then we find 
C2 ~ 0.03. C,2 is small because it is proportional to ^ /rriT^. 

At NNLO, (5=0 once we impose C4 = CIjCq. The condition in Eq. ([u]) 
must still be satisfied. The NNLO good fit condition involves ^3 and C4, 



C3 



2^C +f^V— !^ 



in-^ V87r/^ 

tI, , (Mg\ 
ml^^ \8Trp 



Air 7 
^Mm^ 
An 



Re Li 



27 



12 



(12) 



21n2- (1 



21n2)-^ + of^ 



Since (^4 is multiplied by j'^ /ml, this condition basically fixes the value of ^3. 
At this order C,i may change from its value at NLO. We have chosen to fix C,i 
and (^4 by performing a least square fit to the lower momentum data. 

The -^So phase shift is showp_in Fig. |^. The solid line is the result of 
the Nijmegen phase shift analysig23. The ^Sq phase shift has an expansion in 
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Figure 7: Fit to the ^ Sm phase shift 5 emphasizing the low momentum regim. The solid 
line is the Nijmegen fitEj to the data (for p > 400 Me V values from the VPIE3 phase shift 
analysis were used.). The long dashed, short dashed, and dotted lines are the LO, NLO, and 
NNLO results respectively. AS is the difference between these results and the solid line. 



powers of Q, and we have plotted the LO, NLO and NNLO results. The LO 
phase shift a.t p ^ m^ is off by 49%. At NLO, the error is 14%. At NNLO, 
the error in the ^Sq channel is less than 2% at p ^ m^, and the NNLO result 
gives improved agreement with the data even at p ^ 400 MeV. 

Using M = 939 MeV and m^ — 137 MeV, the parameters for our fit in the 
^Sq channel are: 



NLO 

NNLO 



Ci = 0.2163; 
Ci = 0.0777; 



C2 = 0.0318; 
C2 = 0.0313; 



C3 = 0.1831; C4 = 0.2447; 



Note that Ca ~ Q is larger than (2 ~ Q° because from Eqs. (|l]) and ( [121) 



C3/C2 



7(7A 



NN) 



The parameter (2 is stable because it is fixed by the 



NLO good fit condition. On the other hand, Ci changes by a factor of 2.7 going 
from NLO to NNLO. One expects the value of coupling constants to change 
at each order in the expansion, but a factor of three difference is somewhat 
surprising, ft is also disturbing that ^4 is greater than ^1, since, on thj^ basis 
of the RGE, it is expected that C4 < Ci- At NNLO the RGE for C2 is: 






A 
2/2 



C2 
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Expanding C2 perturbatively results in two equations: 



(13) 






Co,0 + ^^2 ) ^"2 + Co C2,-l 



with solutions 

C2 = (1^0, 



.C2 



Co 



9\ 



^2,-1 — 2—- Co,o + 7775 + C4 Cg 



2/2 



(14) 



The second term in C2,-i has exactly the same form as the leading C2. C2,-i is 
supposed to be a perturbative correction to C2, so C4 ~ Q < Ci ~ Q^j ^ind one 
does not expect this part of C2 _i to be significantly larger than the leading 

C2. 

Some insight into this puzzle can be obtained by comparing the NLO and 
NNLO expressions for tq. Using Eq. (^ we have 



r, 



NLO 




8^ 



1- 



87 272 



NNLO 
1^0 



87r 
M 
8tt 
M 
Stt 
M 



(Ci 



-0.2952), 



SniTr 



C4 + Cl 1 



2TO7r — 27 



A 



NN 



(C4 



0.01375 Cl + 1.35078 C2 + 0.210398) 



47 



(15) 



(16) 



where I/Aatjv — {M g\) / {8tt f'^) . Any reasonable fit for the phase shifts will at 
least approximately reproduce the observed effective range. In Eq. (|l6|), the 
piece of the NLO correction proportional to (i is almost exactly cancelled by 
the NNLO correction. This cancellation occurs because 1 — 2(to^ — 7)/AArjv — 
0.01 in the ^Sq channel. This is simply an unfortunate accident; in the '^Si 
channel, where 7 = 45.7 MeV instead of — 7.88MeV, the coefficient of Ci in 
Eq. ( [iq ) is « 0.4. Since C2 is small due to the NLO good fit condition, this ac- 
cidental cancellation forces ^4 to make up the observed effective range. There- 
fore, ^4 is much larger than anticipated. 

It is important to note that the coupling C2 is not changing nearly as 
drastically at each order. If one were to solve the theory exactly, one would 
find that C2 had a term f Cq , where f represents a short distance contribution 
to the effective range. Ci and (4 can be thought of as the first few terms in 
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an expansion of f. The theory should eventually converge to the exact f but 
it need not reproduce f exactly at low orders in perturbation theory. It is 
reassuring that Cf ^^ = 0.22 and (^^i^^o _^ (^nnlo ^ o.32, indicating that the 
coupling constant C2 is not changing more than one would expect in a theory 
with an expansion parameter ^ 1/3. Note that Q^^^*^ must be small in order 
for c^NNLOj^(^NNLO ^^ ^^^ ^^Q ^^^Yi larger than Cf ^°. At NLO potential pions 
make up ~ 60% of tq, with short distance physics making up the remaining 
^ 40%. At NNLO the situation does not change by very much; potential pions 
give ^ 40%, short distance physics ^ 50% and cross-terms make up the rest. 

Kaplan and SteelecaEj have proposed a fitting procedure in which f is not 
expanded in powers of Q. This amounts to imposing the additional condition 
C4 = 0, so there is no new parameter at NNLO. Only the linear combination 
Ci + C4 appears in C2. However the amplitude depends on C,i and C4 very 
differently because they appear at different orders in the Q expansion. This 
is why we treat Ci and ^4 as seperate parameters. Where it not for the can- 
cellation in r^^^^ noted above, then the difference between the two methods 
would be small, i.e., the size of a N^LO correction. In fact, it is impossible 
to reproduce the observed effective range, if one demands C4 = 0. In the ^Si 
channel their is a logarithmic divergenceQ at order Q, introducing a \ti{^/ K) 
dependence into the coupling C2._i. Since the constant K is undetermined, C4 
cannot be set to zero, so there is a new parameter at NNLO. Fot-Cf this type 
of ln(/i) dependence occurs at order Q^ from soft pion graphsE3 Kaplan and 
Steele have suggested that the failure of their fitting procedure when applied 
to models with effective ranges close to that seen in nature may indicate that 
To is unnaturally large, and that the power counting of the effective theory 
might need to be modified to take, this scale into account. It seems more likely 
that the failure observed in Ref.Ed may just be the consequence of a numerical 
accident in the amplitude at NNLO as shown in Eq. ([iq). This cancellation 
does not occur in the ^Si channel at NNLO, nor is it likely to persist at higher 
orders. For this reason, it seems premature to conclude on the basis of the 
NNLO amplitude that the expansion is failing due to a large rg. 

In Ref. Ej, Rupak and Gautam use a similar fitting procedure to the one 
discussed here. Instead of finding Ci and C4 by fitting to the phase shift, they 
fix these constants by matching onto 

pcoi{5)^-f+'-^{p'+^^) + ..., (17) 

which is similar to the effective range expansion except pcot(^) is expanded 
about p — ij. At NLO Ci is fixed to give sq. At NNLO the same value of Ci is 
used and (^4 is again fixed to reproduce sq. This procedure was applied to data 
from a two- Yukawa toy model and the convergence of the EFT looks similar to 
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that in Fig. M. In our approach wc have not demanded that the exact value of 
rg is reproduced since we know that there wiU be corrections to tq from higher 
orders in the m^/A expansion. 

FinaUy, we would like to comment on the prediction of higher order terms 
in the effective range expansion 

p cot((5) = h -—p^ + V2p'^ + v^p'^ + vap'^ .... 

a I 

Using the NLO expression for pcot(J), Cohen and HanserE3 obtained predic- 
tions for U2,W3 and W4. At NLO, the effective field theory predictions for V2^ 
V3, and W4 disagree with the Vi obtained from a fit to the Nijmegen phase shift 
analysis. The NNLO predictions for the shape parameters are shown in the 
table below. The prediction for tq is not better at NNLO than at NLO. The 
NNLO Vi predictions depend on ^1 and C2. We see that the NNLO correction 
substantially reduces the discrepancy between the effective field theory predic- 
tion and the fit to the Nijmegen phase shift analysis, but the discrepancy is still 
quite large. This gives some evidence that the EFT expansion is converging 
on the true values of the Vi, albeit slowly. Effective field theory predictions 
for the shape parameters have been studicdJjji toy models where one is able to 
go to very high orders in the Q expansion.^ In the toy models, the effective 
field theory did eventually reproduce the shape parameters, but the observed 
convergence is rather slow. 

^■5*0 ■-■ To V2 V3 Vi 

FiS 2.73 fm -0.48 fm^ 3.8 fm^ -17 fm^ 

NLO 2.65 fm -3.3 fm^ 19 fm^ -117fm^ 

NNLO 2.63 fm -1.2 fm'^ 2.9 fm^ -0.7 fm^ 

What can we learn from Fig. |j about the convergence of the KSW expan- 
sion? It is pleasing to see a NNLO calculation reproducing the ^Sq phase shift 
at p ~ 771^ with an accuracy of a few percent, and giving an improved fit to 
the data even for larger momenta. It is important to keep in mind that the 
NNLO calculation of the phase shift is incomplete since there is a possible con- 
tribution from order Qf radiation pion graphs. Many other process involving 
two nucleons can be examined at this order. Once enough processes are cal- 
culated to NNLO, all parameters of the theory appearing at this order can be 
extracted and it will be possible to make predictions with no free parameters. 
The accuracy of these predictions will constitute a serious test of the KSW 
expansion method. 
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Appendix 

Here we give the definitions of the constants appearing in the NNLO expression 
for the amplitude in Eq. (|9|) . After solving the renormalization group equations 
at this order one finds that all quantities in parentheses and curly brackets are 
separately fi independent. The quantities in curly brackets vanish at NNLO 
in the Q expansion. 
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Ci and ^4 are short distance constants of integration of the RGE's in Eq. (y_3|). 
On the other hand, C2 and C3 are sums of constants and variables that appear 
in the amplitude. Note that the order Q radiation pion contribution from order 
Q\ graphs is constant and appears in ^3. AZ)2(/i) is a correction to D^ which 
cancels the ln(/i) dependence from the Q\ radiation pion graphs. 
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